lonization-induced asymmetric self-phase modulation and universal modulational 
instability in gas-filled hollow-core photonic crystal fibers 



Mohammed F. Saleh/ Wonkeun Chang, ^ John C. Travers,^ Phihp St. J. Russell/'^ and Fabio Biancalana^ 

^Max Planck Institute for the Science of Light, Gilnther-Scharowsky str. 1, 91058 Erlangen, Germany 
^Department of Physics, University of Erlangen- Nuremberg, Germany 
(Dated: March 2, 2013) 

We study theoretically the propagation of relatively long pulses with ionizing intensities in a 
hollow-core photonic crystal fiber filled with a Raman-inactive gas. Due to photoionization, previ- 
ously unknown types of asymmetric self-phase modulation and 'universal' modulational instabilities 
existing in both normal and anomalous dispersion regions appear. We also show that it is possible 
to spontaneously generate a plasma-induced continuum of blueshifting solitons, opening up new 
possibilities for pushing supercontinuum generation towards shorter and shorter wavelengths. 
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Introduction — The invention of the photonic crys- 
tal fiber (PCF) has led to a true revolution in the field 
of nonlinear fiber optics [1 . Hollow-core PCFs (HC- 
PCFs) with the so-called kagome-lattice cladding have 
become in recent years a superior host for the investi- 
gation of light-matter interactions between intense ultra- 
short optical pulses and gaseous or liquid media |2l |3] . 
These fibers are characterized by a broadband transmis- 
sion range with low group velocity dispersion (GVD), and 
a high confinement of light in the core. 

In a series of recent groundbreaking experiments, few- 
/iJ femtosecond-scale pulses have been launched into 
argon-filled cm-long HC-PCFs, leading to the unique 
phenomenon of the photoionization-induced soliton self- 
frequency blueshift [4]. These experiments are well de- 
scribed by a unidirectional pulse propagation equation 
(UPPE) based on the full electric field, complemented by 
sophisticated tunneling and multiphoton ionization mod- 
els OE]. However, the UPPE model can be reduced to a 
set of two coupled equations for the envelope of the elec- 
tric field and the ionization fraction [7l|8], and remarkable 
qualitative agreement obtained. 

In this Letter we show theoretically that when the gas 
is excited by relatively long pulses with ionizing inten- 
sities, new kinds of self-phase modulation (SPM) and 
modulational instability (MI) can emerge during prop- 
agation. Moreover, after the initial stage of instability is 
over, a 'shower' of hundreds of solitons, each undergoing 
an ionization-induced self- frequency blueshift, pushes the 
supercontinuum spectrum towards shorter and shorter 
wavelengths. Such a blueshifting plasma-induced contin- 
uum has some similarities with the redshifting Raman- 
induced continuum driven by the Raman self-frequency 
shift in conventional solid-core fibers pHTT]. 

Governing equations — The propagation of light in a 
HC-PCF filled with an ionizable Raman-inactive gas can 
be modeled by the following normalized coupled equa- 
tions [71 E]: 
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where is the electric field envelope, (/) is a quan- 
tity proportional to the number of electrons created 
by the photoionization process, ^ is the normalized 
longitudinal coordinate along the fiber, r is time in 
a reference frame moving with the input pulse group 
velocity, D{idr) is the full GVD operator, a = 
K ((/)T - (/>) [l - IV^Ith/IV^P] (^IV^P) is the ionization- 
induced loss, is a normalization factor, (see also Refs. 
[71 m [12] for extensive details), 6 is the Heaviside step 
function, A|?/^p = — IV^I^^' IV^Ith normalized 
ionization threshold intensity, a and 0t are constants 
that can be determined through the tunneling ionization 
equation rate [3 [8] , and r is a coefficient regulating the 
recombination between free electrons and ions, which is 
a two-particle process and hence depends on 0^. 

Equations ([l][2| yield extraordinary qualitative agree- 
ment with experiments using ultrashort input pulses [3], 
accurately predicting the dynamics of soliton formation 
and the plasma-induced soliton self-frequency blueshift 
[71 [8] . At the core of the model is the assumption that 
tunneling photoionization is dominant in the dynamics, 
and that the ionization rate can be described by a linear 
function corrected by a Heaviside function that takes into 
account the threshold intensity [3l [71 IS! • Although the 
long pulses considered here are expected to drive some 
multiphoton ionization, comparisons of Eqs. (l]2) to 
UPPE simulations using the more general Yudin-Ivanov 
ionization model [13] - which we present at the end of 
this Letter - show excellent agreement. 

Plasma-induced asymmetric self-phase modulation — 
lonization-induced SPM can be studied by simplifying 
Eqs. ([l][2| for the case of small dispersion and long input 
pulse durations to- In this situation (natural for kagome- 
PCFs since the GVD is very small), the nonlinearity ini- 
tially dominates over the GVD, since the 'soliton number' 
N (i.e. the input pulse energy) is large and the sec- 
ond term in Eq. ([T]) can be safely neglected in the very 
initial stage of propagation. In kagome-PCFs, input en- 
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ergies such that N > 100 are reahstically achievable for 
pulses with > 0.5 ps duration. Neglecting also the re- 
combination process, Eq. ^ can be written as (j) (r) = 

{l-exp[-a/:^A|^/;(T')pe(A|^/>(T')P)dr']}. For 

small values of a, Eqs. (Tj2) can be reduced to a single 
integro-differential equation, 



J — oo 



(3) 

where 77 = ck^t- By temporarily ignoring the losses 
(which do not change the qualitative picture that we are 
going to describe, but only saturate the SPM spectrum 
after a certain distance), Eq. (|3| can be solved by sub- 
stituting i/j (^, r) = V (^, r) exp [ (/^nl (^', r)] [12 , where V 
and cp-^i, are the amplitude and the nonlinear phase of the 
pulse. Equating the real and the imaginary parts, and 
performing the integrations results in (^, r) = 7/^ (0, r), 

^NL r) = e - A|F(rOP 6(A|nrOP) 

and A|l/p = |l^p-|^lth- We wih show below that such a 
nonlinear phase can induce a strong spectrally asymmet- 
ric SPM at a large rate compared to the background Kerr 
nonlinearity. For this purpose, we compute the mean 
frequency (1^), and the variance {Any = (n^) - {ny 
[nmi], where (n'') = J n'' d^/ J I'^f dn, = 
■F[fp], and the symbol T represents the Fourier trans- 
form. 

In the case of an input Gaussian pulse, (0, t) = 
exp (— Tq), the above parameters can be determined 
in closed forms. We find 
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2|^l?herf(r) 



(4) 



18 



{d + iVse -3ev^T§ [3erf2 {V2T) 



-2 v^erf (v^r)] + el^Alt^ [erf (r) - 1] 
c [mleTiiT)-V2eTi{V2T)]}, 



(5) 

where erf is the error function, T = T/tq, and — T < r < 
T is the regime within which the pulse intensity exceeds 
the threshold intensity. 

Panels (a,b) in Fig. [l] depict the spatial dependence 
of the mean frequency and the standard deviation, re- 
spectively, for different values of r^, i.e. for different free- 
electron densities generated in the fiber. In Fig. [ija), 
for 77 = 0, which corresponds to the absence of ioniza- 
tion, the mean frequency is always zero during propaga- 
tion due to the well-known symmetric spectral broaden- 
ing due to conventional SPM [12 . As r] increases (i.e. 
when the plasma starts to build up inside the fiber), 
the mean frequency moves linearly towards the blue-side 
of the spectrum due to the ionization-induced phase- 
modulation. This induces a strong spectrally asymmetric 
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FIG. 1. (Color online). Spatial dependence of (a) the mean 
frequency (Q) and (b) the frequency standard deviation of 
a Gaussian pulse exp (— r^/2ro) with tq = 2. The temporal 
position T at which the pulse intensity can initiate photoion- 
ization is assumed to be equal to tq. AQq is the spectral 
width at ^ = 0. 



SPM that is unique to the kind of gas-filled PCF waveg- 
uides studied in this Letter. By looking at Fig. [ijb), 
we find that the ionization-induced SPM broadens the 
spectrum significantly faster than the well-known Kerr- 
induced SPM. The spectral broadening process is limited 
by the unavoidable ionization and fiber losses, similar to 
Kerr SPM. 

Plasma-induced modulational instability — After the 
very initial SPM stage (described above) is over, the in- 
terplay between nonlinear and dispersive effects can lead 
to an instability that modulates the temporal profile of 
the pulse, creating new spectral sidebands referred to 
as modulational instability (MI) in Refs. p!5HT7] ). For 
instance, during the propagation of a continuous wave 
(CW) signal in an anomalously dispersive Kerr medium, 
perturbations in the amplitude of the carrier wave can 
trigger the generation of spectral sidebands, leading to 
the eventual break-up of the CW signal into a train of 
pulses UK [151 nil • MI due to the photoionization non- 
linearity described by Eqs. ( p!p| ) can be investigated by 
using the standard approach presented in Ref. [12, il5j . 
The steady state solutions of Eqs. (lj2), which are time 
independent, are found by inserting i/j = i/jq e*^°^ and 
(j) = Perturbing these steady state solutions and sub- 
stituting back into Eqs. (Tj2), one obtains the linearized 
equations 



drU = 



0, 



(6) 
(7) 



where a and u are perturbations added to the amplitudes 
ijjQ and (/)o, s is the sign of the second-order dispersion 
(+1 = anomalous dispersion, — 1 = normal dispersion). 
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higher-order dispersion coefficients are neglected, F = 
^ {^0 - ml) (V'o' - ml) +2r<j>o, 9 = a{4>T - 0)^0- 

and 
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Assuming a = ai e*'^ + a2 e~*'^ ^ u = uqc 
1^ = — Qr, a set of two homogeneous equations for ai 
and a? can be found, 
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where mn = — k — s ri^/2-|-i/;Q-|-(/i/'o/ {—F + i Wi2 = 
^0 + 5 ^o/ {-F + i m2i = -mi2, and m22 = -Wii - 
2k. This set has a nontrivial solution when 
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which are the complex eigenvalues of the problem. Note 
that in this analysis (j) is kept to be a real positive quantity 
since it represents the number of electrons generated via 
photoionization. 

Perturbations can only be amplified during propaga- 
tion for frequencies which have non-real propagation con- 
stant Hi. The well-known Kerr-induced MI (e.g. in solid 
core optical fibers or in HC-PCFs when excited with in- 
tensities below the ionization threshold for the gas) oc- 
curs only in the anomalous dispersion regime [12]. How- 
ever, the presence of the photoionization process induces 
a very unusual instability that can exist in both nor- 
mal and anomalous dispersion regimes, and for any fre- 
quency. The reason for this 'universal' instability is that 
the ionization contribution to the destabilization of the 
background wave gi^o/ {—F -\-irt) is always complex and 
frequency-dependent, in contrast to when only the Kerr 
effect is present. The spectral dependence of the gain, 
which is defined as 2Im{/^} for different peak powers, is 
shown in Fig. |2]for (a) anomalous and (b) normal disper- 
sion regimes, where the physical powers are normalized 
to the threshold ionization power Pth, i.e., = 1. 

When the normalized input power ipQ < IV^I^h^ we have 
the traditional side- lobes, which exist uniquely in the 
anomalous regime, due to the Kerr-nonlinearity. How- 
ever when ipQ > IV^I^h? photoionization-induced instabil- 
ity causes the side lobes to be unbounded and to have 
slowly-decaying tails. A similar situation occurs in the 
normal regime, where the gain is slightly lower due to 
the absence of the Kerr contribution. For this case there 
are no instabilities below the threshold power since no 
plasma is generated. 

Numerical simulations — To validate the analysis 
based on the simplified model given by Eqs. ( |1|2[ ), and 
also to provide further insight to the ionization driven MI 
dynamics, UPPE simulations were performed [5 . The 
ionization rate was calculated using the Yudin-Ivanov 
model [13 , which accurately describes both quasi-static 
tunneling and multiphoton ionization along with the 
transition between these regimes [131 HI- With this 
model previous experimental results have been repro- 
duced [IS]. 
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FIG. 2. (Color online). MI spectral gain profile versus nor- 
malized input peak power in (a) anomalous and (b) normal 
dispersion regimes with equal magnitude of GVD. 



A Gaussian pulse at 1064 nm with duration 1.18 ps 
and peak power 200 MW was launched into a kagome 
HC-PCF, core diameter 20 /im, filled with argon gas to a 
pressure of 1 bar. This corresponds to anomalous disper- 
sion with (32 ^ -2.8 psVkm. In this fiber Pth ^ 105 MW 
and the pump pulse has a soliton number N :^ 102. The 
peak of this pulse corresponds to a Keldysh parameter of 
7 = 0.785, and for most of the pulse 7 < 5; this is the 
non-adiabatic tunneling regime [13]. 

Fig. [3]^a,b) shows the temporal and spectral evolution 
of the pulse through the fibre. The first stage of propaga- 
tion shows asymmetric spectral broadening towards the 
blue due to ionization-induced SPM as described above 
- see Fig. [sjb) and inset. Immediately after the SPM 
stage, dispersion starts to play a role, and due to the 
combined Kerr and ionization Mis, unbounded side lobes 
are generated and amplified quickly - see the oscillations 
in the time domain in Fig. |3ja). In the third and final 
stage in the propagation, strongly blue-shifted solitons 
are emitted. For comparison. Fig. |3jc) shows propaga- 
tion without ionization. In this case weak side-lobes due 
to Kerr driven MI are visible, but no blue-asymmetry 
is evident. Fig. [sj^d) shows the spectral propagation 
obtained by solving the coupled equations Eqs. (Tp). 
All of the qualitative features of the UPPE simulations 
are reproduced, validating this model and the subsequent 
analysis discussed in this Letter. 

Further insight into the dynamics can be obtained 
from Fig. [4| where the evolution of the cross-frequency- 
resolved optical gating (XFROG) spectrograms of the 
pulse at different positions along the fiber is shown. The 
pulse is initially asymmetrically chirped to high frequen- 
cies in the centre of the pulse, [Fig. |4jb)], due to the 
higher plasma density created at the peak intensities. 
This is the plasma-induced SPM chirp described above. 
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(a) Temporal evolution 
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FIG. 4. (Color online). XFROG spectrograms for the propa- 
gation through an Ar- filled HC-PCF. The simulation param- 
eters and color bars are the same as in Fig. [s] (a) Reference 
pulse, (b) Frequency chirping and initiating MI. (c) Pulse 
break-up. (d) Soliton disintegration into multiple blueshift- 
ing solitons. 



FIG. 3. (Color online). Spectral (a) and temporal (b) in- 
tensity evolution of a long Gaussian pulse propagating in an 
Ar-filled HC-PCF, calculated using a UPPE. (c) The corre- 
sponding spectral intensity evolution in absence of ionization, 
(d) The spectral intensity evolution obtained by solving cou- 
pled equations Eqs. (Tp). 



At the same time two ionization-induced MI sidebands 
appear in the pulse spectrum [Fig. [4jc)], as discussed 
above. MI facilitates the formation of many solitons. In 
less than half a meter of propagation the initial pulse 
disintegrates into a 'shower' of solitary waves, see Figs. 
[4]^d), each undergoing a strong self-frequency blueshift 
induced by the intrapulse photoionization, as described 
in detail in Kefs. fZ',^. 

Conclusions — In this Letter we have investigated 
analytically and numerically the propagation of intense 
and relatively long pulses in an Ar-filled kagome-cladding 
HC-PCF. Several surprising results have emerged, such 
as: strongly asymmetric SPM, which could be studied an- 
alytically with the envelope equations ([l][2|; a universal 
type of MI, existing in both anomalous and normal dis- 
persion for all frequencies and exhibiting long tails in the 
gain spectra; and the final disintegration of the pulse into 
a multitude of blueshift ing solitary waves, which form a 
plasma-induced continuum. All these regimes are real- 
istically accessible in experiments. These theoretical re- 
sults further highlight the stimulating new possibilities 
opened up by accessing the ionization regime in hollow- 
core photonic crystal fibers. 



[1] P. St. J. Russell, Science 299, 358 (2003). 

[2] P. St. J. Russeh, J. Lightwave Technol. 24, 4729 (2006). 

[3] J. C. Travers, W. Chang, J. Nold, N. Y. Joly, and 
P. St. J. Russell, J. Opt. Soc. Am. B 28, All (2011). 

[4] P. Holzer, W. Chang, J. C. Travers, A. Nazarkin, J. Nold, 
N. Y. Joly, M. F. Saleh, F. Biancalana, and P. St. J. Rus- 
seh, Phys. Rev. Lett. 107, 203901 (2011). 

[5] W. Chang, A. Nazarkin, J. C. Travers, J. Nold, P. Holzer, 
N. Y. Joly, and P. St. J. RusseU, Opt. Express 19, 21018 
(2011). 

[6] P. Kinsler, Phys. Rev. A 81, 013819 (2010). 

[7] M. F. Saleh, W. Chang, P. Holzer, A. Nazarkin, J. C. 

Travers, N. Y. Joly, P. St. J. Russell, and F. Biancalana, 

Phys. Rev. Lett. 107, 203902 (2011). 
[8] M. F. Saleh and F. Biancalana, Phys. Rev. A 84, 063838 

(2011). 

[9] M. N. Islam, G. Sucha, I. Bar- Joseph, M. Wegener, J. P. 

Gordon, and D. S. Chemla, Opt. Lett 14, 370 (1989). 
[10] A. S. Gouveia-Neto, M. E. Faldon, and J. R. Taylor, . 
[11] E. M. Dianov, A. B. Grudinin, D. V. Khaidarov, D. V. 

Korobkin, A. M. Prokhorov, and V. N. Serkin, Fiber 

and Integrated Optics 8, 61 (1989). 
[12] G. P. Agrawal, Nonlinear Fiber Optics, 4th ed., San 

Diego, California (Academic Press, 2007). 
[13] G. L. Yudin and M. Y. Ivanov, Phys. Rev. A 64, 013409 

(2001). 

[14] S. C. Pinault and M. J. Potasek, J. Opt. Soc. Am. B 2, 
1318 (1985). 

[15] A. Hasegawa and W. Brinkman, IEEE. J. Quantum Elec- 
tron.. 

[16] K. Tai, A. Hasegawa, and A. Tomita, Phys. Rev. Lett. 
56, 135 (1986). 



[17] K. Tai, A. Tomita, J. L. Jewell, and A. Hasegawa, Appl. 

Phys. Lett. 49, 236 (1986). 
[18] V. M. Gkortsas, S. Bhardwaj, C. J. Lai, K. H. Hong, 

E. L. Falcao-Filho, and F. X. Kartner, Phys. Rev. A 84, 



